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Briefing Purpose 



• Operational perspective: Present challenges that 
space weather poses for performing conjunction 
assessment mission 

• Research perspective: Present NASA Robotic 
CARA research to date on efforts to characterize 
risk from changing space weather predictions on 
conjunction assessment 
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Agenda 



• Background: 

- NASA Robotic CARA 

- Atmospheric Drag 

• Time Offsets: 

- Empirical Evidence 

- Impact on Conjunction Assessment 

• Conclusions and Questions 
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NASA Robotic Conjunction Assessment 

Risk Analysis (CARA) 

• CARA provides support 
to all operational NASA 
robotic missions 

• Supports 67 missions, 
including 

- Earth Science 
Constellation 

- TDRSS 

- Hubble Space 
Telescope 

• As well as a service to 
other agencies 

- NOAA for POES 
satellites 

- USAF for SBSS and 
DMSP satellites 
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The Conjunction Assessment Risk Analysis mission at NASA 
GSFC is to protect NASA robotic assets from threats posed by 
other space objects while operating in the space environment 
through ensuring domain expertise, a robust concept of 
operations . and an operationally-responsive system to meet the 
expanding needs of the mission area 
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NASA Robotic CARA (con't) 



• A conjunction is defined here to 
mean a local minimum in the 
difference between the position 
components of two trajectories - 
the closest point of approach 

• Conjunction Assessment [CAj is 

the process of predicting 
conjunctions 

• Conjunction Assessment Risk 
Analysis [CARAj is the process of 
analyzing the conjunction event 
to determine the associated risk 
to the asset 

• Collision A voidance is the process 
by which the risk associated with 
a conjunction event is mitigated - 
by a maneuver or other action 


JSpOC Orbital Safety Analysts 

• Perform dote approach screening 

• Generate and deliver OCM data to risk assessment team 

• Perform manual ODon secondary object s 

• Request Increased tasking If necessary 


O 


NASA Conjunction Assessment Risk Analysis Team 

• Perform collision risk assessment analysis 

• Quant tfy collision risk 

• Compute Probability of Collision 

• Perform Pc Sensitivity Analysis 

• Perform Pc Forecasting 

• Interface with JSpOC to ensure data product delivery and tracking 
requests 

• Assist mission O/O with avoidance maneuver strategy 

• Provide recommendations to Mission Management and Flight 
Operations teams 


/ \ 


Mission Management 


Mission Operations Team 

• Evaluate conjunction risk 


• Plan and execute avoidance 

relative toother mission risks 


maneuver strategy 

• Make final maneuver or 


• Provide predictive ephe merit 

mitigation decision 


and uncertainty to risk 
assessment team 
• Provide input and 
recommendat ions to mission 
management 
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Conjunction Assessment: 
JSpOC Process and Products 


Now 


Time of Closest 
Approach (TCA) 


SSN obs 


Primary 

Object 

(operating 
object of interest) 

Secondary 

Object 

(debris or 


Batch OD 


.«* 


Batch OD 


K- 


Propagation using 
predicted atmospheric density 
(HAS DM/DC A model with 
NOAA-predicted F10.7, Ap) 


Time 


Predicted 

Conjunction 

Event 


otherwise) J S P 0C Hi S h Accuracy 

Catalogue (ASW) 

• Orbital Conjunction Message (0CM): 

- Includes both object's state vector and position covariance at TCA 

• Allows computation of probability of collision (PJ 

- May receive multiple OCMs over time from additional Space Surveillance 
Network (SSN) tracking 
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Space Weather and Conjunction 
Assessment: A Notional Event 




25 Jan: first identification of possible conjunction on 1 Feb 

27-28 Jan: P c first increases to level of concern before starting to fall (looking safer) 
29 Jan: Alert of a Coronal Mass Ejection (CME) heading for Earth on 31 Jan 
Spacecraft owner/operator (O/O) wants to know if (and how) CME will impact 
conjunction event 

• Does the new space weather prediction make this event safer or riskier? 

• Might performing a maneuver make the conjunction event worse? 
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Space Weather and Conjunction 
Assessment: General Questions 


• Are observed changes in conjunction event 
consistent with space weather changes? 

• How does changing space weather predictions 
affect conjunction assessment predictions? 

- If state vectors and/or covariances impacted, this 
impacts P c 

- How to enfold space weather uncertainty into 
conjunction assessment? How to communicate this to 
O/Os? 

- Assessment of current risk and mitigation strategies: 
too conservative, not conservative enough? 
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Agenda 



• Background: 

- NASA Robotic CARA 
Atmospheric Drag 

• Time Offsets: 

- Empirical Evidence 

- Impact on Conjunction Assessment 

• Conclusions and Questions 
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Background: Space Weather 
Impacts on Satellites 



• Particle bombardment 

- Electrical charging 

- Ionization events 

• Satellite disorientation 


• Communication loss 


Important to 
conjunction - # 
assessment 



a.i. solutions 


Increased atmospheric drag 

- Satellite position 

- Covariance size 

Ionospheric effects 

- Incomplete/inadequate ionospheric correction will impact 
range performance of SSN ground-based radar tracking 

- Could impact OD 


Focus of NASA Robotic CARA space weather 
research relates to atmospheric drag 
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Background: Atmospheric Drag 



• Atmospheric drag magnitude: a dra 9 = \fipv 2 

- n = -2^ is ballistic coefficient 

m 

p is atmospheric density 

v = r'sat 

- Solar cycle and space weather have strong impact on 
neutral atmospheric density 

- Increasing atmospheric drag impacts: 

• Frequency of "Drag Make-Up" maneuvers to stay in control 
box 

• Covariance size 

- Uncertainty in predicted atmospheric drag impacts: 

• Future satellite position predictions (next slide) 
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Background: Atmospheric Drag and 

Predicted Satellite Position 



predicted 

®drcig 



• Satellite will be at a different position if 


n actual predicted 
a drag ^ a drag 


- Uncertainty in predicted atmospheric density not currently 
incorporated into propagation results 

- Uncertainty in ballistic coefficient incorporated in 
covariance 


• Drag acceleration ~counter to satellite velocity 

- Change to drag primarily results in offset in along-track 
position 

- Equivalently can be represented as an early/late offset time 
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MASA ROBOTIC CAR4 



Atmospheric Drag and Covariance Size 


Analytic covariance growth model : 

- Drag case [assuming no uncertainty in p)\ 

34 


Mean anomaly: a 


M 


2 


= f— ) 

a 0 \p o' 


Quadratic growth 
in time 


„ ■ ■ ■ _ .. a 0^0 

• Semi-major axis: &a — — ) 1 

where 

i = d drag ■ v = -ij Bpv 3 = energy dissipation rate 
t = propagation time 
subscript 0 refers to epoch 

Higher drag in the past [during OD) leads to larger 
covariance size in future [at TCA) 


* Extension of Hoots (AAS 11-579) 
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Atmospheric Drag and Covariance Size 

(con't) 

• Effect of larger covariance on P c : 

- At a snapshot of time: Most P c values decrease but some 
P c values (at small misses] increase* 

- As a function of time: Could delay determination of 
conjunction event being assessed as threat/non-threat 
until closer to TCA 

• "Classic” P c time series curve for a miss slowly rises at first but 
rapidly falls-off near TCA as the covariances contract 

• Potential impact to conjunction mitigation timeline 

- Desirable to have clear recommendation before maneuver "go/no- 
go” decision 




* Jenkin (AIAA 2002-1810), Frigm and McKinley (AIAA 2010-7823) 
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Agenda 



• Background: 

- NASA Robotic CARA 

- Atmospheric Drag 

C> Time Offsets: 

- Empirical Evidence 

- Impact on Conjunction Assessment 

• Conclusions and Questions 
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Time and Energy Offsets 
Between Predictions 

• Predicted satellite location at TCA 



(not to scale) 



Can calculate time and energy differences with 
multiple OCMs: 


A f ■ V 

- Time difference: At = , , (rectilinear motion assumption) 


v 


- Energy difference calculated using (osculating) specific 
energy: v 2 \i /i 
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Time and Energy Offsets 
Between Predictions (con't) 

• Changes to satellite energy and time offset sensitive 
to changes in: 

- Propagated atmospheric drag (as predicted atmospheric 
density already incorporated in propagated state vector) 

- Changing satellite position commensurate with position 
covariance 



Early Late 

- Potential ambiguity in interpreting results 

• Attribution to space weather may be impossible but may be able 
to establish plausibility 
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Empirical Evidence for Time Offsets: 

Preliminaries 



• This presentation only examines At 

- Specific energy plots ‘noisy’ and highly sensitive to 
changes in state vector 

• Time offset plots to follow: 

- Time offset plots in narrow altitude band (650 and 750 
km) which includes the NASA A-train 

• Primary (operating) satellites, 6-7 days from TCA 

• Primary (operating) satellites, 2-3 days from TCA 

• Primary and secondary objects, 5-6 days from TCA 

- x-axis: TCA date in decimal days (2012) 

- At measured relative to last OCM received 
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MASA ROBOTIC CAR4 



©<■ 


Time Offset Plot: 
Primary Satellites, 6-7 days from TCA 
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Satellites highly correlated (tend to all be early or late) 

Very small spread in ballistic coefficient values 
Indicative of differences in atmospheric density 
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Time Offset Plot: 
Primary Objects, 2-3 days from TCA 
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Closer to TCA, magnitude of At gets smaller 

Persistence of some of the correlations seen in previous slide 
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MASA ROBOTIC CAR4 



Time Offset Plot: Primary & Secondary 

Objects, 5-6 days from TCA 
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Secondary objects last tracked 5-6 days from TCA 
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Symbol size proportional to ballistic coefficient value 
^ m Larger ballistic coefficient values have larger time offsets 

® a *‘* solutions Ghrist/DeHart/Newman| IMPACT Workshopl Jan 201 3| 21 



MASA ROBOTIC CAR4 



Time Offset Plot: Primary & Secondary 
Objects, 5-6 days from TCA (con't) 
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— . _ Correlation between sign of At changes for primary and secondary objects 
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Empirical Evidence for Time Offsets: 
Conclusions and Lead-in to Impact on CA 


• Strong supporting evidence of changing atmospheric 
drag impact time offsets 

- Time offsets for payloads small but synchronized 

- Objects with larger ballistic coefficients show more severe 
time offsets as they are more strongly influenced by drag 

• What is impact of a time offset on conjunction 
assessment? 

- Qualitative approach: 3 ‘cartoon’ scenarios 

• Both objects early/late by identical amount 

• "Head-on" event 

• "Non-head-on" event with different offset times 

- Quantitative approach 
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Qualitative Impact of Time Offsets: 

Scenario 1 



Both objects early/late by identical amount: 

* 

Original event / Both objects equally late 

/ 

/ 

/ 

/ 

/ 

/ 

/ 



Slide intended to be viewed 
as custom animation (not 
captured in pdf version) 


Results: TCA shifts, but at TCA the conjunction event looks identical 
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Qualitative Impact of Time Offsets: 

Scenario 2 

“Head-on" event: 

Original event 



< 




> 


Both objects have time offsets (not necessarily identical) 


Slide intended to be viewed 
as custom animation (not 
captured in pdf version) 




Results: TCA shifts, but at TCA the conjunction event looks identical 
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Qualitative Impact of Time Offsets: 

Scenario 3 



“Non-head-on" event with different offset times: 


Original event 



Slide intended to be viewed 
as custom animation (not 
captured in pdf version) 


Both objects late by differing 
amounts 



Results: TCA shifts, and conjunction event looks different at new TCA 

© a. i. solutions 


Ghrist/DeHart/Newman| IMPACT Workshopl Jan 201 3| 26 



Quantitative Impact of Time Offsets: 

Theory 


• Miss component changes: related to time offsets 

- Use state vectors from OCM: r p , v p , r s , v s 

- WLOG, shift secondary object by At [net difference) 

- Rectilinear motion assumption: 

= % + v p t jV(t) = v v 


rs(t) = r s +v s (t+At) v ' s {t) = v s 

- New TCA: Find new t [relative to original TCA) 


W 


F , I? 

f rel lJ rel 


= (rj - Vp) ■ (v' -v^) = 0 
= (r s + v s (t + At) - Vp- Vpt ) ■ (v, 

= frel/^rel + &tV s ■ (v s - V p ) + t 

-At (v s ■ (v s - v p )) —At(y s ■ v rel ) 


v) 


V, 


v s - V v 


= 0 


t = 


\v s - V v 


I ^rel I 
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Quantitative Impact of Time Offsets: 

Theory (con't) 






Miss vector at new TCA: 




At(v s ■ Vygl) y _ 

|.v„l ( ”' 



lyel C Vs * r re ; ) I 


Dependent on 
event geometry 


Strategy: vary At to generate family of possible misses 


• Decompose into RIC miss components: 

- Radial: R = f 

- In-track: t= (r x v) x r 

- Cross-track: C = fxv (Normal to orbital plane) 
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Quantitative Impact of Time Offsets: 

Example 1 



I,C zero crossing at ~0.5 sec 
0 d.i. solutions l R l onl y can increase 


SAFE with respect to time offset 
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Miss Component [m] 



Quantitative Impact of Time Offsets: 

Example 2 



There exists a time offset that is extremely dangerous. Might that time offset occur? 
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Quantitative Impact of Time Offsets: 

Results 



• RIC miss components change in a coordinated manner 

- Varies based on conjunction geometry 

• Varying time offset generates family of possible misses 

- Covers how conjunction might evolve over time 

• No predictive power - we don't know the actual time offset 

- Is there a time offset that would result in a potentially dangerous 
situation? 

• Does the radial miss get close to zero? 

• Do all miss components cross zero at around the same time offset? 

• In process of evaluating utility of time offsets for evaluation 
conjunction assessment risk 
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Conclusions and Questions 



• Using time offsets at TCA is a candidate technique to 
quantify how conjunction event might change 

- Time offsets are a physical effect on satellites from changing 
atmospheric drag predictions 

- Some conjunction events at substantially higher risk if there 
is a time offset 

• Critical to give accurate risk assessments in light of 
changing space weather predictions 

• Can we say anything about expected At values for 
changing space weather? 

- Can At be predicted or at least bounded? 

- Can the small At's of primary satellites be used as a 
predictor for the remainder of the catalogue? 
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